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ABSTRACT 


A historical  data  file  based  on  the  near-surface  water  r.asses 
of  the  North  Atlantic  Ocean  is  discussed.  The  most  attractive  feature 
of  the  water  mass  file  is  that  the  characteristics  of  the  input 
bathythermogram  will  objectively  determine  the  proper  deep  r.istory  for 
computation  of  the  surface  to  bottom  sound  speed  profile.  A 
second  feature  is  the  adjustment  of  salinity  to  the  presence  of 
temperature  inversions  (sound  channels)  to  maintain  a stable  water 
column.  Evaluation  of  the  water  mass  file  using  salinity-temperature- 
dcpth  (STD)  data  shows  thcit  is  is  superior  to  the  file  presently  used 
in  the  Integrated  Carrier  TVntisubmarine  VJarfare  Prediction  System  (ICAPS)  . 
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ERRxVTA 


Tables  1,  2,  and  3:  Columns  headings  should  be  changed  from  "RMS 
SD"  to  "Mean  RMS". 
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Sound 


a necessary  input  to  sonar  range  prediction  model; 


Because  synoptic 


sound 


ynthetic  profiles  are  constructed  either  hy  combining  a synoptic  bathy 


thermograph  trace  with  deep  historical  oceanographic  data  (Mendenhall 


Faucher,  et  ai;  Hanssen  and  Tucker)  or  by  historical  data  alone  (Russell 


Podeszwa) . Each  of  the  several  techniques  available  relies  on  a 


different  method  of  generating  the  surface  to  bottom  sound 


They  agree,  however,  in  that  they  provide  a 


ingle  seasonal  profile  for 


each  region;  with  each  region  having  fixed  boundaries.  Tinfortunately 


real-world  oceanographic  features  are  constrained  on 


by  bathymetric 


boundaries,  and  their  position  may  vary  rcipidly  as  in  the  case  of  a 


Gulf  .Stream  meander  or  cyclonic  eddy.  Thus  historic  files  based  upon  a 


single  regional  history  frequently  provide  misleading  sonar  range 


oceanographic  data  file  — based  on  water  ma 


in  which  the 


bathythermograph 


computer  program  uses  the  characteristics  of  the  inpu 


trace  to  automatically  select  one  of  several  possible  histories.  The 


file  was  designed  to  he  incorporated  into  the  Integrated  Carrier  Anti 


submarine  Vlarfare  Prediction  System  (ICAPS)  developed  under  the  cognizance 


U.S.  Naval  Oceanographic  Office  (NAVOCEANO) 


Two  i)remises  were  made  while  developing  the  new  file:  (1)  that  near 


surface  water  masses  can  be  uniquely  identified  by  thermohaline  characteristics 


fuid  (2)  the  l•.hermal  characteristics  of  noi ghbc)ring  v/att.-r  masses  arc 
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ficiently  different  so  as  to  permit  reliable  identification  from  an 


expendable  bathythermograph  (XbT)  trace  alone.  After  identification  o: 


the  applicable  deep  history,  temperature  values  of  the  input  trac 


inarged  with  deep  temperatures  using  an  equation  of  the  form 


where  T.  and  Th.  are,  respectively,  estimated  and  historical  temperature 


at  depth  i,  K a weighing  factor,  and  AT  the  difference  between  temperature 


at  the  bottom  of  the  XBT  trace  and  interpolated  historical  temperature 


at  the  same  depth.  The  weighing  factor,  developed  from  inperical 


solution  for  a set  of  historical  data,  is  determined  as  a function  of  the 


depth  increment  between  points 


Because  few  guidelines  for  water  mass  identification  are  available  in 


it  was  decided  that  th 


method  of  determining  water  mass  characteristics  within  a given  area  was 


to  look  at  original  oceanographic  data.  Two  world-wide  data  files  were 


approximately  491K  observations  compiled  by  the  National  Oceanographic 


Data  Center  (NODC)  provided  temperature  and  salinity  data  at  each  of  32 


expendable  bathythermograph  (XBT)  file  of  approximately  218K  observations 


comoilcd  from  three  sources  (NAVOCEANO,  NODC,  iuid  the  Fleet  Numerical 
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Weather  Ceetor)  provided  toiaperature  data  at  each  flecture  point  over  tlie 
depth  range  of  the  iiistrinnent  (as  deep  as  760  m)  . The  following  procedure 
was  used  to  determine  water  mass  characteristics  in  the  near-surface 
layer  (0-400  m) : 

A.  The  classical  literature  was  searched  for  applicable  descriptive 
papers.  For  example,  the  northern  edge  of  the  Gulf  Stream  is  frequently 
delineated  by  the  15°C  isotherm  at  200  ra. 

B.  The  ocean  .station  data  file  was  used  to  provide  annual  composite 
statistical  data  (mean,  standard  deviation,  number  of  observations)  at 
each  standard  depth  using  all  available  data  within  the  area  of  interest. 

A plot  of  the  distribution  of  temperature  versus  salinity,  plotted  at  both 
200  and  400  m,  provided  insight  as  to  the  number  of  water  masses  present 
and  thermohaline  variability  within  each  water  mass.  Figure  1 shows  a 
plot  of  temperature  versus  salinity  at  200  m in  the  rectangle  45  to  50“n, 

40°  to  50°W  — on  area  where  the  cold  I..ibrador  Current  meets  with  the  warmer 
North  Atlantic  Drift.  The  presence  of  water  masses  with  specific  thermo- 
haline characteristics  are  clearly  recognizable,  and  tentative  water  mass 
classification  has  been  made.  The  200  m level  was  found  to  be  an  ideal 
depth  for  classification  in  that  it  is  generally  well  below  the  level  of 
both  diurnal  and  seasonal  changes  while  being  within  the  depth  range  of 
both  XBT  and  AX3T  (airborne  deployed  XBT)  probes.  The  XBT  file  provided 
statistical  data  and  histograms  for  temperature  and  temperature  gradients 

at  preselected  depths  to  supplement  the  ocean  station  data  when  necessary'. 

C.  Flecture  points  in  the  temperature  versus  salinity  (T-S)  plot 
shown  in  Figure  1 clearly  defined  water  mass  criteria  in  areas  where 
different  water  masses  exist  in  close  proximity.  Considerable  temperature 
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vnriataiiity  also  occurs  in  areas  occupied  by  a single  water  rr.as.s,  probably 
a result  of  dynamic  events  such  as  uo'velling.  iv'here  vari.ibility  of  this 
nature  was  observed,  two  classifications  ('warn'  ojid  'cold')  were  made 
to  provide  a better  merge  between  XliT  trace  and  history. 

D.  Temperature  filters  at  200  in  were  d'veloped  to  distinguish  adjacent 
wator  masses  based  on  information  provided  in  tne  px'ovious  step.  Where 
adjacent  water  masses  had  similar  toniperature  range  at  the  test  depth, 
differentiation  was  made  by  examination  of  the  temperature  gradient 
between  the  200  and  300-m  levels.  For  example,  both  the  Oulf  Stream 

and  the  Sargasso  Sea  are  cliaracteri^ed  by  a temperature  range  of  from 
15“  to  25“C  at  200  m.  A near- Lsotherm.al  layer  of  18“C  v;ater  is  found  to 
extend  from  the  bottom  of  the  seasonal  thermocline  to  depths  exceeding 
300  ni  in  the  Sargasso  Seia.  No  such  layer  exists  in  the  Gulf  Stream. 
Examination  of  Sargasso  water  in  a region  well  removt'd  from  the  Gulf  Stream 
{30°-35°N;  60°-65°N)  showed  that  95  percent  of  the  observations  had  a 
temperature  gradient  between  0.0“/100  m and  -1.6“C/100  :n.  Thus  a gradient 
of  -1.6“c/100  m at  the  200-300  m level  is  used  in  the  region  of  the  Gulf 
Strea.m  to  differentiate  Sargasso  Water  from  Stream  Water. 

E.  Mean  seasonal  temperature  and  salinity  values  were  then  determined 
for  each  depth  and  water  mass  (figure  2).  Where  the  data  are  not  sufficiently 
daetj  temperature  and  salinity  were  extrapolated  to  the  bottom  by  comparison 
v;ith  neighboring  profiles.  Inconsistencies  in  the  data  — such  as  temperature 
inversion  at  depths  below  200  m — were  examined  to  determine  if  they  are 

a result  of  statistical  processing,  data  distribution,  or  bad  data. 
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Temperature  and  salinity  at 
standard  depths  in  Slope  VJater 
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F.  A quality  control  check  was  made  by  plotting  the  sranonal  data  { 

on  a single  plot  of  tempecwature  versus  salinity  (figure  3)  . Inconsiatancies 
in  the  data  are  immediately  apparent;  temperature  errors  by  a vertical 

spilce,  salinity  errors  by  a horizontal  spike,  and  depth  errors  by  a skewed 

r 

spike.  Where  data  were  obviously  incorrect  the  rjlot  was  smoothed  to  con- 
form with  surrounding  data. 

i 

EVALUATION  j 

The  file  was  evaluated  by  comparing  the  merged  ijrofile  generated  by  both  | 

the  now  water  mass  file  and  the  old  file  — which  is  based  on  a single  j 

seasonal  profile  for  each  5-degree  rectangle  --  v/ith  oceanographic  data.  ! 

I 

Test  data  typical  of  each  \i;ater  mass  within  the  test  areas  v,-e.re  selected  j 

I 

from  salinity-temperature-depth  (STD)  observations  on  file  at  cither  ilODC 

or  the  Coast  Guard  Oceanographic  Unit.  Six  observations  per  season,  divided 

equally  among  water  masses,  were  selected  for  each  area.  The  uppermost  j 

I 

^ portion  (0-400  m)  of  the  STD  cast  was  treated  as  an  XBT  and  the  temperature  j 

trace  extended  to  1,500  m by  merging  with  both  old  and  new  history  files.  } 

J 

Salinity  was  estimated  ai\d  sound  speed  computed  for  all  depths  between  the  j 

surface  and  1,500  m.  In  the  surface  layer,  estimated  salinity  and  sound  j 

speed  values  were  compared  with  observed  values  for  each  depth  on  tlie  ; 

simulated  XBT  trace.  In  the  deep  layer  (400-1,500  m)  estimated  temperature, 
salinity,  and  sound  speed  values  were  compared  v/ith  observed  values  at  6 ' 

depths:  500,  600,  800,  1000,  1200,  and  1500  n. 

The  first  tost  was  designed  to  test  the  premise  that  quality  controlled 
data  from  a large  area  — in  this  case,  a 5 x 10-degree  rectangle  — would 
compare  favoraibly  v;ith  the  smaller  area  without  quality  controlled  data  as 
used  to  compile  the  old  history.  Should  this  promise  prove  correct,  then 
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considerable  reduction  could  be  made  in  the  file  size.  The  area  including 
Ocean  Weather  Station  EGiO  (44'‘N,48°W)  was  selected  because  considerable 
STD  data  were  available  from  an  area  of  relatively  little  ocecuiograpaic 
variability.  The  results  of  the  test  at  OWS  ECHO  are  given  in  Table  1 . The 
new  v;ater  mass  file  provided  slightly  better  results  in  both  the  svirface 
and  deep  layers. 

The  second  test  was  designed  to  document  the  ability  of  the  new  file 
to  dif Eerentiate  between  mater  masses,  thereby  providing  a merged  profile 
superior  to  that  produced  by  the  old  file.  An  area  of  high  oceanographic 
varicibility  seaward  of  the  Virginia  Capes  (VAC7VPES)  was  selected  because  of 
the  presence  of  3 water  masses:  Slope  Water,  Gulf  Stream  Water,  and  Sargasso 
Water.  Results  of  this  test  (Table  2)  show  that  the  water  mass  file  estimate.s 
salinity  significantly  better  in  the  surface  layer  with  a corresponding 
increase  in  the  accuracy  of  sound  speed  computations.  The  deep  data  again 
are  slightly  better  v/hen  estimated  by  the  water  mass  file  than  with  the 
old  file. 

The  final  test  evaluated  the  ^Jbility  of  the  water  mass  file  ho  adjxist 
salinity  values  in  a near-surface  temperature  inversion  (sound  channel) . 
Persistence  of  sound  channels  for  months  at  a time  show  that  they  are 
stable  oceanographic  features.  However,  use  of  unadjusted  historical 
salinity  values  causes  an  apparent  unst.ible  water  column.  Thus  historical 
salinities  must  be  reduced  by  the  method  given  in  the  appendix  if  they  are 
to  bo  realistic. 

Salinity  adjustment  was  evaluated  in  Slope  Water  in  the  VACAPKS  area 
whore  v/ell-def ined  inversions  occur  from  April  through  October.  Salinity 
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was  estimated  in  an  inversion  using  the  water  mass  hLstori’  first  with  and 
then  without  the  adjustment  routine  and  compared  with  observed  values  from 
20  STD  drops.  Results  of  the  evaluation  — given  in  Table  3 — indicate 
that  adjusted  salinity  values  are  clearly  superior  to  unadjusted  values  in 
temperature  inversions. 

The  new  water  mass  file  offers  svibstantial  advantages  over  conventional 
oceanographic  data  files.  The  major  advantage  of  the  riew  file  is  the 
capability  of  the  system  to  automatically  dlscem  in  which  of  several 
possible  ’.•/ater  masses  an  input  observation  was  taken,  thus  providing  for 
variable  water  mass  boundaries.  An  evaluation  of  the  v/ater  mass  selection 
capability  showed  that  significantly  better  results  wore  realized  by  its 
utilization.  A second  feature  of  the  new  file  is  salinity  adjustment  within 
a temperature  inversion.  A test  of  this  capability  showed  that  is  provided 
a realistic  estimate  of  salinity  in  an  inversion.  Finally,  the  quality 
control  procedure  to  which  the  water  mass  file  is  subjected  permits  the  use 
of  large  areas  without  loss  of  accuracy. 

Several  recommendations  are  warranted  based  on  the  evaluation  of  the 
water  mass  file: 

A.  The  water  mass  file  should  be  expanded  to  include  t'ne  Indian  and 
North  Pacific  Oceans. 

B.  That  the  V'/ater  mass  file  replace  the  5-degree  history  file  presently 
used  in  ICAPS. 

C.  That  the  water  mass  file  should  be  adopted  as  the  navy -wide 
oceanographic  history  file  for  on-scene  applications. 
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D.  Salinity  in  tho  proposed  file  is  dstiimiined  as  a function  of  water 
mass,  season,  and  depth.  In  future  water  laass  files  it  is  racoiraiended  that 
s.alinity  be  determined  as  a function  of  water  mass,  temperature,  depth, 
and  season  to  take  advantaije  of  the  unique  tliermohaline  relationships  w)iich 
by  definition,  identify  any  given  water  mass  using  a technique  similar  to 
that  used  in  the  western  North  Atlantic  Oce;in  by  f’isher. 

E.  Because  the  effect  of  averaging  masks  real  featiires,  it  is 
recommended  that  .in  ancillary  file  containincj  digitized  X3T  traces  typical 
of  each  water  mass  on  monthly  basis  be  constructed  to  provide  users  with 
infonnation  vital  to  realistic  planning  of  future  fleet  operations. 
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SALINITY  ADJUSTMENT 


Oceanographic  stability  is  a prime  requisite  for  persistence  of  near- 
surface temperature  inversions,  for  unstable  conditions  v/ouli  destroy  an 
inversion  through  mixing  in  a relatively  short  period.  In  physical 
oceanography,  stability  is  quantified  by  the  change  in  density  with  respect 
to  change  in  depth;  stable  conditions  being  denoted  by  a positive  gradient. 
Because  density  is  a function  of  both  temperature  and  salinity,  it  follows 
a priori  that  a salinity  inversion  must  coincide  with  a temperature  inversion 
if  stability  is  to  be  maintained.  Historical  salinities,  by  definition, 
represent  mean  conditions  and  thus  cannot  cope  with  an  anomalous  condition 
such  as  an  inversion.  This  appendix  describes  a method  of  adjusting 
salinity  as  estimated  from  an  historical  file  to  provide  a stable  water 
column.  In  order  to  allow  for  minor  instabilities  frequently  in  Arctic 
waters,  the  correction  is  only  applied  where  a temperature  inversion  exceeds 
0,25°C. 

The  equation  used  to  adjust  historical  salinity  ’.';as  derived  from  stepwise 
regression  of  density  as  a function  of  salinity  (30  to  40  “/oo)  a.iid  constant 
temperature  {10“c) 

p = -1.26584xlo''^+7.72412xlo"^S+4.22003xlO~®s'^  (1) 


Differentiation  of  equation  (1)  to  give  change  of  density  v;ith  respect 
to  c;hange  in  salinity  yields,  after  rearrangement,  addition  of  a correction 
term  to  assure  stcibility  within  the  inversion,  and  conver.sion  of  density  to 
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the  more  conventional  sigma-t  yields 


As  ■ = 

0. 7724+1.6880x10“ 'S  3 
o 


(2) 


i'Anere  Aat  : tlie  difference  between  Cf^  at  adjacent  points, 

: the  corrective  term,  and 
Sq  : original  historical  salinity  at  point  i 

The  initial  step  in  applying  equation  (2)  is  the  computation  of  sigma-t 
as  a function  of  depth  and  temperature  as  input  from  the  XBT  and  interpolated 
historical  salinity.  The  XBT  is  scanned  from  bottom  to  top  and  salinity 
values  adjusted  for  all  points  within  temperature  inversions  that  are  more 
than  0.25°C  less  than  the  temperature  ma.ximum  at  the  lower  boundary  of  the 
inversion.  The  corrective  term  is  increased  by  0.01  sigma-t  units 
per  depth  increment  to  assure  stability  v/ithin  the  inversion.  7idjusted 
historical  salinity  (Sj^)  is  now  computed  using  the  equation 


Si  = + ASi 


(3) 
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